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and the Heavy-light Spectrum from NRQCD 

A. Ali Khan a * 

a Physics Department, The Ohio State University, Columbus, OH 43210, USA. 

The present status of lattice calculations of the B spectrum and /s, using NRQCD for the b quark, is discussed. 

1. INTRODUCTION 



B physics is a subject of active research. The 
theoretical and experimental understanding of 
the B meson and b baryon spectrum is just be- 
ginning. Weak matrix elements of B mesons, 
e.g. fs, fs s , Bb, and Bb s are being stud- 
ied to determine fundamental parameters of the 
Standard Model. This review summarizes the 
progress made in calculating the B spectrum and 
decay constants using the nonrelativistic QCD 
(NRQCD) [| approach on the lattice. 

The advantage of NRQCD is that the rest mass 
term is removed from the Lagrangian. Hence, 
large O(Ma) effects pose no problem, and one 
can simulate b quarks directly on the lattice. 
Alternate approaches have been reviewed by 
T. Onogi Q| at this conference. 

In heavy- light mesons, NRQCD is equivalent 
to a 1/M expansion. The rationale is that if 
the heavy quark is nonrelativistic (p = Mv), and 
the light quark is relativistic with a momentum 
p ~ Aqcd, then momentum conservation in the 
meson rest frame gives 

Mv ~ Aqcd- (1) 
For B mesons, one has: 

Aqcd/M~ 0.1. (2) 

It is thus reasonable to include relativistic cor- 
rections in an expansion in powers of Aqcd/M . 
The lowest terms in the expansion are as follows. 
The 1/M corrections to the static NRQCD La- 
grangian density L — D t are 
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rif =2 staggered lattices from HEMCGC 
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16 3 x 48 



Table 1 

NRQCD calculations of the B spectrum and de- 
cay constants. For references see text. 
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*Talk presented by A. Ali Khan at LATTICE '97, Edin- 
burgh, July 1997. 



In the following, actions which only include H^ 1 ' 
will be referred to as 1/M actions, and those 
which include H^ 2 \ as 1/M 2 actions. Note that 
7j( 2 ) also includes the first relativistic correction 
to the kinetic energy which is formally 0(1/Af 3 ), 
but is expected to give a contribution of a similar 
size as the spin-dependent interactions 0{1/M 2 ). 
An overview of the simulations covered is given 
in Table [j]. For lack of space I have not included 
a recent calculation on coarse lattices using im- 
proved glue and light fermions gj. 

In this review, I will make comparisons to ex- 
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periment and to other lattice results where pos- 
sible. To judge the reliability of the predictions, 
I also evaluate systematic effects like the effect 
of truncating the 1/M expansion, quenching er- 
rors, and the dependence of the decay constant 
on lattice spacing. 

2. SPECTRUM 

Using nonrelativistic b quarks and clover light 
quarks on the lattice, it is possible to reproduce 
the presently known general features of the B 
spectrum. An overview of the the meson spec- 
trum at a lattice spacing a -1 ~ 2 GeV is pre- 
sented in Fig. [I]. The plot shows results from 
the most comprehensive calculations on quenched 
(1/M 2 action, GLOK §j) and dynamical (1/M 
action, SGO ||) configurations. The B s — Bd 
splitting agrees with experiment, also the B(2S) 
agrees well with the first experimental candi- 
date §. 

For the P wave states, two kinds of experimen- 
tal signals have been found. A £?Wtt resonance 
has been established which is expected to be a su- 
perposition of various P states. The states with 
light quark angular momentum ji = 1/2 (Bq and 
B{) and ji = 3/2 (Bl' and B^) are expected 
to form doublets. The splittings within these 
doublets are given by the coupling of the heavy 
quark spin to ji . The second experimental signal 
is a narrow Btttt resonance which has probably 
31 = 3/2 §. 

Lattice NRQCD predicts a B^ - Bq splitting 
of the order of 200 MeV (an effect of - 4a). The 
B\ and B* 1 cannot yet be resolved since each of 
the two operators used in the simulation for the 
j = 1 states project only on a superposition of 
these states, called B{ in Fig. [I]. For rif — 2 lat- 
tices, there is presently only a P state signal in the 
channel that (similar to the quenched data point 
on the left of it) corresponds to the 1 P\ state in 
heavy-heavy mesons. Given the limitations that 
apply to both the experimental and lattice results 
on P states, my conclusion is that there is quali- 
tative agreement. 

A study of the effect of the 1/M 2 terms in 
the action on some spectral quantities is ongo- 
ing Q . A low statistics comparison between 1 /M 




Figure 1. The B meson spectrum from lat- 
tice NRQCD. Circles denote quenched results, 
squares, results from n/ = 2 lattices. Dashed 
lines denote experimental error bounds as given 
by the Particle Data Book, dotted lines the cen- 
tral values of first experimental measurements of 
these states. 



and 1/M 2 actions on quenched configurations at 
(3 = 6.0 d,[§, suggests that there is no significant 
difference. Assuming that this is true, a compar- 
ison of the spectrum on quenched and the nj = 2 
lattices from HEMCGC, as shown in Fig. [j], fur- 
ther suggests that quenching errors are small too. 

Lastly, it is important to study the dependence 
of the results on the lattice spacing a. Prelim- 
inary results of a comparison between quenched 
lattices at (5 = 5.7 and 6.0 show no significant 
scaling violations ||. 

2.1. The Hyperfine Splitting 

The B* — B splitting is of particular interest 
since it is expected to be sensitive to various ef- 
fects such as discretization errors, quenching and 
the tuning of the parameters in the heavy quark 
action. The results shown in Fig. ^ are ~ 20 MeV 
too low. The question arises whether the results 
improve for smaller a. Fig. || shows results us- 
ing Wilson and clover light quarks at various lat- 
tice spacings, with a determined from m p . With 
present errors, I conclude that there is reason- 
able scaling for clover light quark data at j3 — 5.7 
and 6.0, as well as for Wilson light quark data 
at (3 = 5.8 and 6.0. Also, there is no obvious 
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Figure 2. NRQCD results for the B hyperfine 
splitting. Crosses use Wilson light quarks, circles 
tadpole-improved clover light quarks. The solid 
and dotted lines denote the experimental value 
with its error bounds. 

difference between clover and Wilson results. 

If the B hyperfine splitting depends on the 
heavy quark wavefunction at the origin, one 
would expect quenching effects to be significant. 
Comparing results on quenched and dynamical 
(rif — 2) configurations (Refs. HQ], does not 
show any significant effect. One possibility is 
that the physical picture of a nonrelativistic wave- 
function is incorrect in B mesons. However, the 
sea quark mass in the HEMCGC configurations 
might be too large, so modern dynamical simula- 
tions with high statistics are needed to settle this 
issue. Another quantity whose effect on the hy- 
perfine splitting needs to be studied are the coeffi- 
cients of the spin-dependent terms in the action. 
A perturbative calculation by H. Trottier |jj is 
under way. 

2.2. Comparison to other lattice calcula- 
tions 

Fig. H shows data from NRQCD and other lat- 
tice methods along with a preliminary experimen- 
tal result for a P state excited B s meson. In some 
of the calculations only an approximate value for 
the strange quark mass is used |lC|Jll|| . Also, the 
calculations were done at slightly different (3 val- 
ues between 5.9 and 6.1. (The lattice spacing was 
set using either m p with clover quarks or yfd. 
These give roughly consistent values for a fixed 




Figure 3. Compilation of lattice results for ex- 
cited B s mesons. Circles are NRQCD results, di- 
amonds pi and triangles |11| show static heavy 
quarks. The square shows tadpole-improved 
clover heavy quarks 0|. Results for the Bq and 
the f?2 are grouped together wih dotted circles. 
The meaning of the dashed lines is the same as 
in Fig. [p. At the bottom of the figure, the or- 
bital quantum numbers of the respective states 
are shown. 



(3 in this (3 range.) Because of these limitations, 
a check of the expected 1/M corrections (O(50) 
MeV) to the static limit cannot be made with 
these data. Lastly, I would like to mention that 
first results for D and F wave states are reported 
in Ref. @. 

2.3. b Baryons 

Lattice results for b baryons with non-static b 
quarks are summarized in Fig. |J. For the Af, the 
quenched lattice calculations using NRQCD @, 
tree- level clover Jl3| , and plain Wilson heavy 
quarks agree within errors with experiment and 
with each other. The A^ from Ref. || is signif- 
icantly heavier than the experimental value and 
heavier than the quenched NRQCD result. This 
needs further study. 

For the E& and baryons there are only pre- 
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Figure 4. Lattice results on baryons contain- 
ing one b and two light quarks. Only lattice 
calculations with moving (non-static) b quarks 
have been considered. Circles denote NRQCD 
b quarks j|, diamonds show clover (csw = 1) 
heavy quarks |13|, and the triangle, Wilson heavy 
quarks Q. The square is from a SGO NRQCD 



calculation with n 



f 



The meaning of the 



lines is the same as in Fig. |[ 



liminary experimental results available fl5|| , but 
they are in agreement with the available lattice 
results (see Fig. ||). Ref. also has results for 
strange b baryons. 

3. DECAY CONSTANTS 

At tree level, the nonrelativistic current, cor- 
rect through 0(l/M 2 ), is given by: 



A) = 97o75 Q 

- ]^7o75 (lJ->) Q 



(5) 
(6) 



+ ^97075 [f> 2 +gtB-2iaE)Q, (7) 

where a = 707. 

The question is whether the truncation of the 
l/M expansion at this order is sensible. At 
present, the effects of including both l/M and 
l/M 2 corrections have been studied in two cal- 
culations: Ref. []l6| with Wilson light quarks and 
Ref. Q with clover light quarks. First I exam- 
ine the effect of the l/M 2 terms in the action 
on the matrix elements of the zeroth order cur- 



Figure 5. Ratio of the unrcnormalizcd current 
corrections (left, 0{1/M), right, 0(1/M 2 )) to the 
uncorrected current. Crosses use Wilson light 
quarks at (3 = 5.8 circles, clover light quarks 
at j3 = 6.0 (GLOK). For the x axis, an arbitrary 
numbering is used. 

rent (Eq. (|)): Comparing the l/M and l/M 2 
actions with clover light quarks at = 6.0, one 
finds a difference of ~ 7% [§,0- However, this is 
only a ~ 2a effect as the results from JlT]] are low 
statistics. Similarly, the Hiroshima group finds 
the effect of including the l/M 2 corrections to be 
< 0.5% @. 

Next the effect of the l/M 2 current corrections 
is studied. Fig. || shows ratios of the l/M and 
1 /M 2 current corrections to the zeroth order cur- 
rent (Eq. ||) at the mass of the B meson. The 
ratio of the 0(1/ M) operator is ~ 10%, and of 
the 0(1/M 2 ) is ~ 2%. After including renormal- 
ization factors in the ratios, the l/M correction 
is still ~ 10% [Q, whereas the renormalization 
constant for the 1 /M 2 current corrections has not 
been calculated yet. 

The renormalized matrix elements are obtained 
from a matching calculation to full QCD. For 
NRQCD heavy quarks, two such 1-loop calcu- 
lations for operators through l/M exist. In 
Ref. jl9) the full operator mixing at 0(a) has 
been calculated for clover light quarks. This cal- 
culation shows that there is a significant mixing 
with an 0(ak.QCDa) current correction. Ref. pq| 
presents a calculation for Wilson light quarks 
which takes only the diagrams with vanishing ex- 
ternal momenta into account. These give only the 
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Figure 6. Comparison of fs s from NRQCD, us- 
ing configurations with nf = and nt = 2 
Quenched results using a 1 /M action |t7| (circle) 
and a 1/M 2 action (J] (square) are shown. The 
lattice spacing has been fixed from m p (a -1 ~ 2 
GeV). For better comparison, the 1/M 2 current 
corrections are not included in the result with the 
1/M 2 action. 



contributions of Eq. (g) and Eq. (||) on the renor- 
malization of the zeroth order current Eq. (p|) . 
The neglected terms are however significant Jl8[ . 
Also note that for Wilson light quarks the full 
mixing calculation at one loop contains an in- 
frared divergence unless the O(aAQCDa) correc- 
tion is neglected. 

In Table |2|, a summary of fs is given; NRQCD 
data includes renormalization constants. I give 
two results for the Hiroshima group: (i) from 
Ref. using the renormalization constant de- 
termined by them, and (ii) using the full per- 
turbative Z's (but neglecting any 0(o.A.qcd<i) 
terms) at 0(1/ 'M) |19| on their raw data from 
Ref. [Q. To obtain a renormalized value for fg 
from the raw data provided by Draper et al. [22|, 
I also use the renormalization constant from [19|. 
However, note that these Z's are for a slightly 
different action than used in Refs. Jr|^2| ■ 

The first error bar in the NRQCD results in Ta- 
ble H is statistical (including fitting uncertainties 
where applicable). The second error estimates 
the uncertainty in the determination of a -1 . For 
clover light quarks, the upper error bound comes 
from determining a -1 from f v ; the lower bound 



from either ^fa Wl or from the statistical error 
in m p j|. Ref. 0] used an a" 1 of 2 GeV with 
an error of 0.2 GeV. For Wilson light quarks, the 
upper error bound comes from the error in a -1 
from m p , the lower error bound from the deter- 
mination of a -1 from y/a. Here, the difference be- 
tween a^ 1 from m p and y/a at (3 — 6.0 is ~ 20%, 
whereas for the clover action they are roughly 
consistent. The third error bar is an estimate of 
higher order contributions in perturbation theory, 
obtained by using ap Q evaluated at q* = 1/a 
and at q* — -k /a in the Z's, and, for calculations 
including only 0(1/ M) corrections it includes an 
estimate of the 1/M 2 contributions. 

Table |2| also lists results on Jb from relativis- 
tic heavy quark actions. For more details on this 
data see the review by T. Onogi It is en- 
couraging to note that the result with NRQCD 
heavy and clover light quarks [Q agrees well with 
those results with relativistic heavy quarks where 
extrapolations to a — > have been done [p5|-p7|. 
The results with NRQCD heavy and Wilson light 
quarks |J|,|§ in Table | are high (> 200 MeV). 
This is mainly due to the choice of the lattice 
scale; using a -1 from yfa instead of Wilson m p 
yields a result close to §|,§5|-f27|] . Wilson light 
quarks give rise to considerable discretization ef- 
fects O(Aqcd), and if they are used in conjunc- 
tion with NRQCD heavy quarks, these effects 
cannot be removed by extrapolating to a — ► 
(which is in principle possible e.g. in p7[). 

Since fg from NRQCD cannot be extrapolated 
to a — > 0, it is important to understand the dis- 
cretization effects. There exist only preliminary 
results for Jb s as reported by J. Hein et al. ||. 
They find slight scaling violations: fg s decreases 
by ~ 2a between (3 = 5.7 and 6.0. To make a clear 
statement about discretization effects, at least 
one other calculation at smaller a is required. 

Another important question is how quench- 
ing affects the heavy-light decay constant. In 
Fig. ||, fs s from two quenched calculations fll7j[i[ 
is compared with a result using nf = 2 lat- 
tices jy. To accommodate the fact that the 
strange quark masses in the two quenched calcu- 
lations have been determined by either setting K 
or K* mesons to their physical value, the system- 
atic error from various methods to fix the strange 
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Table 2 

Summary of results on fg from moving (non-static) heavy quarks. For the results marked with *, I have 
included the renormalization constants from Ref. p9|| . "t.i. clover" stands for tadpole-improved clover. 
The results from Refs. [n6p2|3,H,E6| and the systematic errors from |25| are preliminary. 
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quark mass have been included in the error bars. 
Existing data in Fig. ^| show no significant devi- 
ation between the quenched and the dynamical 
results. 
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